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Abstract: A model for the dd — > amr reaction, based on two parallel AW — > 
dn reactions, is extended to incorporate a complete set of input amplitudes. While 
low-energy cross sections are underestimated, the rich structure observed in the 
a-particle momentum distributions for 0.8 < T d < 1.9 GeV (the ABC effect) is 
extraordinarily well reproduced. In addition, a recent measurement of deuteron 
analyzing powers agrees quite well with our predictions, both in frequency and 
magnitude of the oscillations. 
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1 Introduction 



On many occasions, interesting physics has been obtained from detailed studies of 
momentum distributions from nuclear reactions. Excitement once arose from the 
unexpected width of /3-decay spectra, which led Pauli to postulate the existence 
of the neutrino. In most cases, however, the interest has focused on details of 
the distributions, such as an enhancement or a suppression at a certain missing 
mass, or skewness, or other irregularities. These deviations from phase space 
mark the discoveries of new particles and resonances, increased propensity for a 
particular partial wave, or some other peculiarity. Regardless of the problem at 
hand, further investigations often proved fruitful and well worth the effort. A 
particular interesting problem of this kind is the striking phenomenon observed 
in several experiments on the pn — > dX Jj], ^, pd — > 3 HeX j|, 0, [5], and 
dd — > aX H |7| reactions. In all three cases there is a peak just above the two- 
pion threshold with a missing mass Mx ~300-320 MeV/c 2 . Since its position 
and width vary with kinematical conditions, it could not be a new particle or 
resonance. The only alternative is some sort of kinematical enhancement in 
double pion production. The experiments reveal that the effect only occurs in the 
isospin Ix = channels of the two first reactions; the dd — > aX reaction leads to a 
pure Ix = state, which is probably the reason for the more pronounced peaks in 
this case. In association with the peaks at low tttt mass there are also broad bumps 
at maximal missing masses. These are most clearly seen in the pn — > dX and 
dd — > aX reactions, which could be due to more symmetric kinematics. These 
characteristics are collectively known as the ABC effect after the discoverers @]. 
The suggestion of the original authors that the peaks are due to a large tttt s-wave 
scattering length was soon ruled out by separate experiments || and theoretical 
calculations [[J. In addition the central bump (unknown at the time) cannot be 
explained in this way, since it would mean that parallel and anti-parallel pions are 
simultaneously favored by the irn interaction, which is very unlikely. The mass 
of the central bump is also changing monotonically with beam energy. It is thus 
impossible to explain the ABC enhancement in terms of some peculiarity of the 
tttt interaction and one has to include effects resulting from the presence of other 
particles in the reactions and consider the structure of the single pion production 
mechanism itself. In all three reactions there are indications that the ABC effect 
disappears at low beam energies [pi], [11], [12], suggesting that the subprocesses 



in any explanation should be strongly energy and angular dependent in order to 
reproduce the quite different dynamics observed in the different energy regions. 
For the pn — *> dX reaction, on which theoretical investigations have mainly fo- 



cused, the most promising approach is the A A- model [14 . Here two independent 



pion productions are achieved through the excitation of both incident nucleons 
into A-isobars via meson exchange, with the deuteron then being formed through 
a final state interaction of the recoil nucleons arising from the A decay. By using 
the energy and angular dependences of A-excitation and decay, the authors could 
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reproduce the shape quite well, but underestimate the cross section at larger an- 



gles ||. A double- nucleon exchange model, also proposed in the seventies JL5 
and using separate diagrams for the three different peaks, showed a much poorer 
angular distribution || and the authors were incapable of obtaining an absolute 
normalization. In his survey of all ABC experiments and models, Barry [|16j tried 
to improve the original AA-model and to extend it to pd — > 3 HeX but with little 
success; his model has, in particular, problems in reproducing the spectra away 
from the ABC peaks. Better results were achieved in a chiral bag model [l7|, tak- 
ing into account 80 diagrams to fourth order in the ttNN coupling constant, but 
this model still failed in the details. A model based on excitation of the W*(1440) 
(Roper) resonance has recently been proposed for the lower energies The 
calculations are in good agreement with data [TIJ , but the contribution from A A 
excitations is negligible. The situation at low energies is thus improving and 
might form a basis for a renewed investigation of the ABC effect in the pn — > dX 
reaction. 

Despite these efforts, there is as yet no completely satisfactory implementation 
of a model for the ABC enhancement in either the pn —>■ dX or pd — > 3 HeX 
reactions. On the other hand, the dd —>■ aX reaction has not been subjected to 
any theoretical investigation because of its perceived complexity. In a preliminary 
study [[0| we were able to reproduce very well the angular dependence measured 
at Td = 1.25 GeV using a version of a A A model where independent single 
pion productions arise from two AW — > dn reactions occurring in parallel. In 
this work, only the dominant AW — > dn input amplitude was kept and we now 
lift this restriction by including a complete set of amplitudes. This makes it 
possible to reproduce all the major features observed in the dd — > aX reaction 
throughout the A region, including also the recent measurement of deuteron 
analyzing powers @. The main idea behind the model is explained in Sec. The 
general description is then given a precise formulation in Sec. |3|, to be followed 
by some remarks about the numerical calculations in Sec. |j. Our predictions are 
compared with data in Sec. |^ and conclusions and outlook are given in Sec. ^. 
Some formulae for the AW — > dn reaction, in particular the relations between the 
partial wave and spin amplitudes, are collected in an appendix. 



2 Model for dd — ► cxttw 

Our model for the dd — > aim reaction is based on a simple, semi-classical picture 
and illustrated as a Feynman diagram in Fig. [|. The two-pion production in the 
dd collision is viewed as two free, parallel, and independent AW — > dn reactions 
taking place between separate pairs of nucleons from the two deuterons. The 
a-particle is then formed by fusing the two final deuterons. 

An intuitive motivation for this model can be obtained by looking at the 
possible configurations of final particles. If the deuteron and a-particle binding 
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energies are neglected, the local NN — > c?7r cm frames obviously coincide with the 
overall cm (CM) so that the two deuterons have the same energy, though not nec- 
essarily the same directions. However, because of the strong p-wave dependence 
of the NN —>■ dn reaction, the deuterons and pions are preferentially emitted at 
small angles to the beam direction. When deuterons are emitted parallel to each 
other, they are easily bound together into an a-particle. Due to the kinematics, 
the associated produced pions then have small relative momenta and hence small 
invariant mass. This corresponds to the situation at the ABC peaks. In the 
second case of anti-parallel deuterons, the pions go back-to-back, their invariant 
mass is maximal and this is the reason for the central bump. However, it then 
needs a lot of Fermi momentum in the initial deuterons and/or final a-particle 
in order to get the deuterons to stick together so that, despite the kinematical 
enhancement of this configuration, it is suppressed compared to the parallel case. 

Both the ABC effect and the central bump are simultaneously explained in 
our model by the strong p-wave dominance in the NN — > dn reaction, with the 
former being helped by favorable kinematics. 



3 Formal description 

The considerations of the previous section are immediately transformed into the 
Feynman diagram of Fig. [l], where the momenta of the particles are defined in 
the CM. In the following paragraphs the different parts of the corresponding 
matrix element are established using a description in terms of non-relativistic 
wave functions. The derivation of phase space formulae is also given. 



3.1 Relativistic phase space 



air 1 7r 



The phase space calculations for the dd 
tivistically starting from the general expression pO 



reaction are performed rela- 



where -Ej,p; and ujj,kj are the energies and momenta of the initial and final 
particles, respectively. The a-particle is detected in the laboratory frame and 
thus |vin| =p/E and d 3 k a /uj a = [k^ dVladka/uj^i^, so that 



f ^^\M\ 2 S\ Pl+P2 -h-k 2 - k a ). 



(3.2) 

The remaining integrand in eq. ( |3.2j ) is relativistically invariant and the re- 
sultant integral is most easily evaluated in the irn rest frame, where the pions 
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go back to back. Denoting with an asterisk quantities evaluated in this frame, a 
straightforward calculation gives 

r c] 3 k* r\ 3 k* k* r 

z= ^^\m^5\k + K)KMx-u{-^) = — <m*\Mr- (3.3) 

The missing mass Mx in the dd — > aX reaction is, in our model, the effective 
mass of the pion-pion system m n7r . In terms of this and k* = |(k* — k^), the 
differential cross section is 

( dV \ 1 k 2 aUh k* f.„l 



dVt a dk a ) l h 32(2tt) 5 m d M x poJa\ab 



JdQ*- ^ \M*\ 2 , (3.4) 



lab "■"V""; "W-'AF ^alab " ^ spin,pol. 

where the averaging over the initial spins and summing over the final is now 
included. 

The description has so far concentrated on 7r + 7r~ production. Simple isospin 
arguments suggest that 7r°7r° production should be exactly half that of the charged 
pion but the narrowness of the ABC peak and the significant pion mass differences 
make it necessary to evaluate the kinematics separately in the two cases. 

3.1.1 Transformation to the nn rest frame 

In order to implement the phase space considerations of the previous section, 
the matrix element has to be expressed in the nn rest frame so that a Lorentz 
transformation is needed between the CM and Tin systems. Since in general the 
a-particle emerges at an arbitrary angle #cm with respect to the beam direction, 
this involves a rotation. The righthanded CM coordinate system is oriented with 
its z-axis in the beam direction (p), the x-axis in the plane spanned by p and the 
a-particle momentum (k a ), and the y-axis parallel to p x k a . A new coordinate 
system is introduced, by rotating about the y-axis, such that the new z-axis 
coincides with the a-particle direction. A Lorentz boost is then applied in the 
a-particle direction so that, in the new system, the pions emerge back to back. 
Explicitly this gives 

u u* 

(3.5) 
(3.6) 
(3.7) 
(3.8) 

where xyz refer to the CM, and x'y'z' to the nn systems and 7^ = (u>i + U2)/m n7r 
is the Lorentz boost factor between the two. The transverse and longitudinal 
components of the pion momentum difference become 



1, 


- ^2) 




\(K 


-k%) 


= k* x , cos 6>cm + 1-kK, sin 9 C m 




-kl) 


k* 


5<*S 




= -k* x , sin 6'cm + 1-kK' cos ^cm 
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f jk^ — k|j = [(A;:,cos^cM + 7A*'Sin^ C M) 2 + (A;;,) 2 ] 5 (3-9) 
\{k{-kl) = —k* x , sin^cM, +1-kK ! cos #cm (3.10) 

where b indicates directions perpendicular to the beam. 

3.2 Vertex parametrizations 
3.2.1 Split-up vertices 

Since large Fermi momenta are not required in our model, we shall only retain 
the dominant S'-state parts of the nuclear wave functions. In a non-relativistic 
pole model the d:pn vertices may be parametrized as: 



M d% = • e,j (2nf 2 ^( q 2 + c&M) (3-H) 

where y?(q) is the deuteron S-state wave function in momentum space, = 
(qf, ?f/7) the relative momentum of the nucleons, boosted in the z-direction to 
take the relativistic motion into account, and 6j the deuteron polarization vector. 
The parameter a d = \JmB d , where m is the nucleon mass and B d the deuteron 
binding energy. 

The dd\a vertex is parametrized in a similar manner by 

M a = • e^j (27T) 3 / 2 2 v / 2^(q 2 + a 2 )^(q a ), (3.12) 

where e and e' are the deuteron polarization vectors, ^(q a ) is the a:dd S'-state 
wave function in momentum space, q a = qi — q 2 — |(ki — k 2 ) is the relative 
momentum of the deuterons, and a 2 = m d B a , with B a being the a-particle 
binding energy. Due to the large mass of the a-particle and the relatively low 
pion momenta, there is no need to consider a Lorentz contraction factor in this 
case. 



3.2.2 Pion production vertices 

The matrix element for the free NN — > dn reaction is expressed in terms of the 
partial wave amplitudes given in the appendix; 

M nt = K tT <T-Q t , (3.13) 

where Ki is the sum over all initial spin-singlet amplitudes and cr • Qj is the sum 
over all triplet amplitudes. Note that, since we are now using NN — > dir rather 
than nd — > NN, the deuteron polarization vectors should be charge-conjugated 
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compared to those in the appendix. The =F signs for the triplet terms, referring 
to the upper and lower vertices of Fig. p] arise from the inversion of the direction 
of the incident nucleon momenta in the two cases. The values of the amplitudes 
were calculated from the SAID database plj using the normalization discussed 



in the appendix. As seen in Fig. |2], the nine amplitudes retained in the appendix 
reproduce the SAID predictions for the ir + d — > pp differential cross section to 
within a few percent. For comparison the cross section calculated purely with 
the l D 2 p amplitude, renormalized to fill the forward cross section, is also shown. 
This single amplitude, which was used in our preliminary study [|l9j, shows the 



manifest pion p-wave dominance which is the main origin for the ABC enhance- 
ment in our model. 

In addition to considering far more input amplitudes, in the present analysis 
we retain also the explicit energy dependence prescribed by the database (see 
Sec. in the earlier work we extracted the threshold k factor and then assumed 
a constant reduced amplitude. The energy dependence of the 1 D 2 p amplitude is 
hence different in the two analyses. 

3.3 Matrix element 

The matrix element deduced from the Feynman diagram of Fig. [I] involves inte- 
gration over two Fermi momenta 

where M tot is composed of the non-relativistic reduced vertex functions and the 
propagators stripped of their spin structure. In an order suggested by the dia- 
gram, 

M tot = ]T M Q P d P d ,Tr(M d2 P 2 c M^ 1 P 4 M dl P3 C M 7r2 P 1 ), (3.15) 

int. spins 

where Pj are the propagators, M x the matrix elements representing the vertices 
and the sum is over the spins and polarizations of the internal particles. Moreover, 
c denotes charge conjugation. 

3.3.1 Approximation of propagators 

The denominators of the nucleon propagators are factorized to make the pole 
structure explicit: 

i 2m 

Pi 4 = 7 '■ tt t — 5 5 : r~, tt (3.16) 

(920,10 +Po ~ «e) [920,io - («d + <E,i ± P • q2,i)/Po + 

— i 2m 

2 ' 3 (920,10 -Po + ie) [920,10 + (a 2 d + qii =F P • QvVpo - ie \ 
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Since the Fermi momenta (q«) are generally quite small compared to the 
incident deuteron momenta, only the poles closest to the origin at gio,2o — + 
Qi 2 -F P • qi.,2) yield significant contributions. The integrations over g 2 o and qio 
are then done by contour integration in the lower half plane, yielding 



d ?2o dgip ~(2m) 4 -m 2 

(2vr) 2 2 4 3 1 (2 Po ) 2 (q? + a3)(q3 + a 2 ) ~ 7 2 (q? + a })(c£ + ' 



where the last step assumes that po = E d ~ 27m. 
The deuteron propagators are 



d,d' 



d) 
(3.18) 



(3.19) 



where Q d ^ = [p ± (?io - <?2o) - ^1,2, ±(qi - q2) - ki )2 ]. Since q 2 + a 2 ~ 
— — m d + Qi' — m d)i the product of the deuteron propagators and the dd:a 
vertex reduces to 



M a P d P d , = I -j=e ■ e'j (2n)*V2^ c 



1 



+ 



1 



C&-m 2 d + ie. Ql -m 2 d + ie / 
The denominators in this expression can be expanded as 

Qld> - m d = [ E d - Wi >2 ± (gio - Q20)? - [±(qi - q2) - ki, 2 - m 2 ? 



(3.20) 



~ =F2m d 

In the last step, the relation 

1 



• (qi + qa) + ]:{oj x - u 2 ) 



Po 



p • (qi + q2) - q? + q| — v d ■ (qi + qa) 



(3.21) 



(3.22) 



has been used, while quadratic terms were neglected. Using these linearized 
forms, the principal value terms cancel in the sum of propagators in eq. (|3.20j ) to 
leave only a 5-function term; 



■1 



M a P d P d , = _ e • e' (2tt) * y/fr^ipl 61 q{ + q z 2 



mdVd 



L0 1 - UJ 2 

2v d 



(3.23) 



Assuming that Fermi momentum effects may be neglected in the spin cou- 
plings, which is a good approximation in the ABC peak regions, the matrix 
element is written as the product 



M = -i— W/C, 

v d 



(3.24) 



where the spin kernel /C and dimensionless scalar form factor W are given, re- 
spectively, by 



K = —^e ■ e'Tr 



mt. spins 



° ' 62 Vi - * ■ Qi) — < ~7^~~ ) (^2 + o- ■ Q 2 ) 



y/2 J V V2 

(3.25) 

and 



3.4 Summing over spins and polarizations 



(3.26) 



The spin kernel of eq. (|3.25|) can be simplified somewhat by expanding the product 



and taking the trace. This results in the spin-decoupled form 

-1 



A(e l ■ e 2 ) - zB • (e x Xe 2 )-e r C- e 2 j , (3.27) 
where the spin-0, -1, and -2 amplitudes are 



A = ^(e-eO^i^ + lQi-Qa) (3.28) 

pol. 

B = 5>-e')(^iQ 2 + ^ 2 Qi) (3.29) 

pol. 

t = £(e-eO(QiQ 2 + Q 2 Qi-fQi-Q 2 ) (3.30) 

pol. 

and the sums are over the polarizations of the internal deuterons. 
The unpolarized matrix element squared is then 

£ |/C| 2 = |(3|^| 2 + 2|B| 2 + p| 2 ), (3.31) 

ext. pol. 

where \C\ 2 = an d the sum is over the polarizations of the external 

deuterons. The deuteron vector and tensor analyzing powers A y and A yy are 
obtained from the sums 

£ KvPjkKl, (3.32) 
where the are the spin projection operators in a Cartesian basis [22||. Since the 



dominant l D 2 p amplitude by itself gives only a contribution to the A amplitude, 
this would lead to A y = A yy = 0. The analyzing powers are therefore very 
sensitive probes of the importance of non-dominant NN —> dn input amplitudes. 

To avoid the tedious algebra associated with nine partial wave amplitudes 
occurring bilinearly in each amplitude, the above sums and contractions were in 
practice carried out numerically. 
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3.5 Form factor in configuration space 



The form factor of eq. (|3.26| ) is most easily evaluated by transforming to config- 
uration space; 



1 f 



x exp 



■-(kx-ka) -x c 



exp 



-- — (^i - lo 2 ){zi + z 2 ) 
4v d 



(3.33) 



where x Q = (b, (zi— z 2 )/2). 

Since the wave functions have been taken to be spherically symmetric, the 
angular integration over the b variable can be performed explicitly. Furthermore, 
only the even parts of the exponentials will contribute, so that 



^yj- j 1 / Z\ -^2 

bdbdz 1 dz 2 $ d (b,jzi)$ d (-b,jz 2 )yUh 



y^h^J ™—V- 5 ( ~x, -»v ~, -ay-, 2 

x J (f |k*-k*|&) cos[|(^ - k\)(z x - z 2 )] cosf^^i - w 2 )(*i + za)][3.34) 

3.6 Reduction of cross sections due to flux damping 

The fluxes of initial and final particles are reduced because of multiple scattering. 
Any detailed microscopic estimate of such effects would rest upon assumptions 
that are difficult to justify quantitatively. For example, in any Glauber-like ap- 
proximation, correlation effects would be very large for the light nuclei that we 
are considering. We therefore limit our ambition, attempting a crude estimate 
through an overall scale factor. This should at least give a hint of the magnitude 
and direction of the damping. Since our model is based upon two NN — > dn 
reactions, it is natural to estimate the reduction in flux through the factor 



V 



2[a tot (pn) + a tot (pp)} 




(3.35) 



where the last factor is squared since there are two emerging pions. 

The assumption behind this formula is that in the initial state the flux reduc- 
tion due to multiple scattering is the same in all dd reaction. A similar argument 
is valid in the final state, assuming the a-particle to be a dd aggregate. The val- 
ues of total cross sections needed were taken from p3| (NN), p4| (dd), p5| (nd), 
and 



TTCt 



Interpolation between different data sets was often needed and 
the values used are collected in Table [1], together with the calculated reduction 
V. The latter is compared to the factor needed to scale our model calculations 
to fit the data, as shown in the figures. While our simplistic estimate goes in the 
right direction, it does not account for all the deviations from the data. There is 
however considerable uncertainty in the pion distortion at the lower energies. 
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4 Numerical calculation 



The evaluation of cross sections and polarization observables was done in two 
steps. The form factor of eq. ( |3.34j ) was firstly obtained using standard Gauss- 
Legendre routines with 50 points in each of the three dimensions and tabulated 
in steps of the parameters fcaiab, f*, and cos6**. In the second step the spin sums 
and contractions were done at each of these points, followed by the integration 
over phase space, including the tabulated form factor. 



4.1 Wave functions 



The iS-state Paris [27] wave function is used to describe the deuterons, whereas 



the a-particle is represented by the dd:ot cluster function of Forest et al. [^8 
It is shown in the latter work that the nucleon-nucleon distributions in the a- 
particle and the deuteron are very similar for short distances (r < 2 fm), with a 
constant scale factor of 4.7 between them. Since pion absorption on deuterons 
occurs mainly when the nucleons are close together, we need the number N a of 
pairs of such "small" deuterons in the a-particle. Assuming the neutron-proton 
distribution inside the two deuterons of the a-particle to be independent of each 
other, we therefore normalize the wave function to N a = 2.3. 

Since the shape of the a:dd function is similar to a Gaussian, exp(— /3 2 r 2 ), it 
is easy to test our dependence on it by varying the parameter (3. The result of 
this test (using the constant amplitudes of next section) at Td = 1250 MeV and 
9 a = 0.3° is shown in Fig. ^. While the ABC peaks remain largely unaffected, 
there is large sensitivity in the central region. Calculations at 11.0° and also 
at 787 MeV show less variation. Hence this uncertainty is associated with high 
missing masses. 



4.2 Partial wave amplitudes 

The values of the partial wave amplitudes can be extracted from the SAID 



database |21] in several ways. The first and simplest is to use constant am- 
plitudes evaluated at the kinematics corresponding to the edges of the spectra, 
i.e. close to the ABC peaks, where the two pions have identical momenta. The 
single amplitude used in our previous work was extracted in this way Jl9| . 



Energy dependent input — where the amplitudes are calculated separately 
for each point of integration — is an attempt to take the Fermi motion in the 
a-particle into account. There is no unambiguous way to do this and so we have 
adopted two slightly different methods. One based on the previous technique, but 
calculated at each point separately, will be called the El prescription. The second 
approach (E2) uses the same pion momenta but now assumes the deuterons to 
have half the a-particle momentum. The two clearly coincide at the edges of 
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the missing mass spectra, but have different characteristics in the central region. 
Their relative merits will be discussed in the next section. 



5 Results 

The predictions of our model are compared with the Saclay unpolarized cross 
section data of Ref. || in Figs [| and [5[ Shown are the raw results from the E2 
prescription and also those obtained by smearing El and E2 over a Gaussian 
experimental resolution in the a-particle momentum with a relative standard 
deviation a{k a )/k a = 0.5% ||. The calculations are scaled to fit the data with 
factors listed in Table [I]. This scaling is independent of the choice of amplitudes. 
The agreement with data is remarkably good, both in angular distribution and 
energy dependence. The differences found between the scale factors and the flux 
damping estimates could be attributed mainly to the crudeness of the estimate 
and the uncertainties in its input data. 

The slight discrepancy at the central bump for small angles could be due 
to the production of three pions or to our simplified treatment of the Fermi 
momenta. A broad cr-meson [23| would also enhance the cross section here. 



An estimate of the possible contribution from production of //-mesons could be 
obtained from a measurement at 1.95 GeV |[29|| . The cross section is however 
so low (~ 1 /ib/(srGeV/c)) that it should be invisible in all the plots, except 
possibly for 1.94 GeV where one might see a hint of a peak. As illustrated in 
Fig. H the central region is also sensitive to the precise form assumed for the 
a-particle wave function. Inclusion of -D-states in the wave functions would give 
effects for large missing mass as well since it corresponds to large Fermi momenta. 

Apart from a small change in the ABC peaks, the main difference between 
El and E2 is found in the central region. The reason for this is the different 
kinematics used in the two cases — El tends to overstate the 7rd energy in the 
central bump and thus attains the maximal AW — > dn cross section already 
at 940 MeV, while E2 has its maximum at 1250 MeV. This accounts for the 
El overestimating (compared to E2) at low missing masses and underestimating 
at high maximal missing mass. At 1250 MeV and 11.0° the missing mass is 
not high enough to be sensitive to the difference between the two assumptions. 



Unfortunately the range of data in the SAID database [^TJ does not extend to the 
energies needed for El at 1412 MeV and for E2 at 1938 MeV. The close agreement 
between E2 and the calculation with constant amplitudes at 1412 MeV justifies 
the use of constant amplitudes at the higher energies in Fig. |5|. 

The underestimate at 787 MeV is probably an indication of a more severe 
problem within our double-A model. There are two measurements at even lower 



energies, one close to threshold at T<j = 570 MeV [[13] and another at 650 MeV [[T2]]. 
In both cases the data are structureless, and consistent with pure phase space, 
with total cross sections of cr(570) ~ 43 nb and <r(650) ~ 600 nb, respectively. 



11 



In our model a pronounced ABC-like structure should be clearly visible at both 
energies. However, our predicted total cross sections are a factor of about 20 too 
low and it therefore seems that there is some other mechanism which is mainly 
responsible for the smooth behavior near threshold. This extra contribution 
might account already for our underestimate at 787 MeV. 

The deuteron analyzing powers measured by the SPESIII spectrometer at 
Saturne j7| are seen in Fig. ^ together with our model predictions integrated 
over the spectrometer acceptance. The slope of vector analyzing power A' y in the 
forward direction is calculated according to the experimental average procedure 
described in Ref. 0. The differences between the El and E2 calculations are 
clearly seen for A' , while they are rather marginal for A yy . This greater sensitivity 
of A' y could be due to its dependence upon the relative phase between singlet 
and triplet amplitudes and also the importance of the strongly energy- dependent 
and dominant 1 D 2 p amplitude. The El predictions are shifted by an arbitrary 
—2 rad -1 in the figure to show the close resemblance to the data structures. Any 
small extra terms in the amplitude could cause such a shift. For both analyzing 
powers, there are again discrepancies in the central region, presumably for the 
same reasons as for the unpolarized cross sections. 

All calculations were done using the C500 solution of but, in order to 
check the stability of the result, they were repeated for the SP96 solution as 
well. The two sets of results are very close and the small differences found are 
completely overshadowed by the larger uncertainties inherent in the procedures 
for extracting amplitudes. 



6 Conclusions 

For the first time a quantitative model has been proposed for the ABC effect in 
the dd — > aX reaction. Despite its simplicity, and without benefitting from any 
free parameters, it is able to reproduce all the main features of the a-particle 
momentum spectrum observed throughout the double-A region. It is hence clear 
that, at least for this reaction, the ABC effect is indeed a kinematical enhance- 
ment in the independent production of two p-wave pions when these emerge with 
parallel momenta. 

The analyzing power predictions are equally impressive and reproduce quite 
nicely the frequency and strength of the oscillations for both A yy and A' y . Since 
these quantities are sensitive to the non-dominant input amplitudes, and in par- 
ticular their relative phases, this more thorough examination of the details of 
our model adds further strength to the conclusion that it supplies the correct 
dynamics in the double-A region. 

The low energy data present, however, a more delicate situation to interpret. 
We do not think that it is possible to force the present model to reproduce pure 
phase space at low energies. Since the p-wave dominance of NN — ► dir extends 
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almost down to threshold, the s-wave pions needed to furnish isotropy could 
not be produced via that subprocess. The analogous contribution of isovector 
nucleon-nucleon pairs, in for example pp — > ppir , is known to be much smaller 
than isoscalar pairs both near threshold and throughout the A-peak |3(J. There 
must therefore be another production mechanism present, which is dominant at 
low energies but which assumes lesser importance at higher energies where the 
ABC is seen clearly. In order to improve our understanding of the dynamics in 
this kinematical region, we suggest that further experiments should be done at 
deuteron beam energies Td < 0.8 GeV. In addition to inclusive polarized and 
unpolarized cross sections, exclusive reactions are also needed so that angular 
distributions of pairs of the final particles can be constructed. Such data already 
exist for pd — > 3 He7r + 7r~ [ITT]. Since all such quantities can be calculated within 



our model they provide extra confirmations and tests of it. 

Given its quantitative success it is natural to ask if our model could be ex- 
tended to include other reactions as well. A simple extension would be to calcu- 
late the double photon production dd — ► 077 via two np — > dj reactions. This 
could be an important background in the measurement of the charge-symmetry- 
breaking dd — ► air reaction, observed at Saturne in the double- A region [ j3l| . 
We are currently investigating this possibility. 
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A Partial wave amplitudes and Clebsch-Gordan 
coefficients expressed by spin amplitudes 



We find it simpler to estimate the double-A contribution to the two-pion pro- 
duction using a spin-amplitude description of the n + d — > pp reaction rather than 
employing partial waves. We therefore give here expressions for the lower partial 
waves in terms of the Pauli spin matrices (cr), the deuteron polarization vectors 
(e), and the proton and pion momenta (p and k). The partial waves are denoted 
by 2Spp+1 {Lp P )jl n , where S pp and L pp are the spin and orbital angular momentum 
of the proton-proton system in the initial state, J is the total angular momen- 
tum, and l n is the pion angular momentum in the final state. The operators 
Oi corresponding to particular partial wave transitions are given below in the 
normalization /dfi£e Tr|C £ | 2 = 2 • 4vr(2J; + 1). 

A.l Singlet amplitudes 




(A.l) 



(A.2) 



O (W) 




A.2 



Triplet amplitudes 




(A.5) 



O ( 3 P 2 d) 




)} (A.7) 



(A.6) 




+ 



(p k)cr • (kxe) + (cr-k)p- (kxe)]} (A.8) 
^{[5(p-k) 2 -l]e-(pxcr) 

10(p • k)(p • e)k • (pxcr) - 2(k • e)k • (px<r)} (A.9) 
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A. 3 Clebsch-Gordan coefficients 



In the cr-e expressions for Clebsch-Gordan coefficients, the non-relativistic two- 
component spinor 77 and its charge conjugate rj c = —ia 2 rj* are introduced. 



1A 



1 1 

-777.1 —mo 

2 2 



-mi 



1 I 

2 

1 MX 1 
— 777.0 1A -7771 

2 1 2 

(oo| W/i 2 > 



-1 



;(P1) . AM) 



(1A| I//1I//2) <-> e^t • f^=e (w) xe^ 2) 



(A.10) 
(A.ll) 
(A.12) 
(A.13) 
(A.14) 



A. 4 Relation to the SAID formalism 

The cr-e formalism uses the normalization 

where ■Js ir d is the cm energy. The matrix element is then of the form 



M. — ( pp 



where .Aj is the complex amplitude for a particular partial wave. 
In the formalism of the SAID database []2I] 



1=1 



(A.15) 



(A.16) 



(A.17) 



where the helicity amplitudes, Hi, are linear combinations of partial wave am- 
plitudes T. The relation between the cr-e and SAID partial wave amplitudes is 
then 



(A.18) 



777 y pk 

since all the spin-momenta structures and relative phases are incorporated in the 
cr-e formulae. 
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Figure 1: Feynman diagram for the dd — > aX reaction showing the momenta in 
the overall cm system 
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Figure 2: Differential cross section for the ir + d pp reaction at T n = 181 MeV, 
calculated with nine amplitudes (dashed line) compared to the differential cross 
section obtained from the SAID database (solid line). The cross sections for 1 D 2 p 
(uppermost dotted line) and the singlet (dot-dashed line) and triplet (lowest 
dotted line) amplitudes are also given. 
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Figure 3: Comparison of cross sections calculated with different a-particle wave 
functions. Solid lines correspond to the Gaussian wave function with the pa- 
rameter (3 (defined in the text) having the values (from bottom to top): 0.55, 
0.64, 0.67, 0.70, 0.75, and 1.0 fm _1 . For comparison the calculation with the 
Forest et al. 



28] wave function (dashed line), is repeated. The data are from 



Ref. 
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V (estimate) 


1.50 


1.38 


0.77 


0.46 


0.48 


0.78 


D(fit) 


2.2 


0.90 


0.72 


0.73 


0.80 


1.0 



Table 1: Distortion factor D calculated according to eq. ( ]3.35|) for different dd 
energies. The total cross section data are from Refs E31 E3L EH |2H] . The estimated 



reduction factor is compared to the scale factor needed to fit our predictions to 
data. 
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Figure 4: Angular distribution dd — > aX at Td = 1.25 GeV. The smeared El 
(dashed line) and E2 (solid line) calculations are compared to the data of Ref. |J. 
For E2 the raw calculations (dotted line) is also given. Our predictions are fitted 
to the forward (low k a ) ABC peak with a scale factor 0.71, cf. Sec. |3.6| . 
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Figure 5: Predictions of the energy dependence of dd — > aX at 6 a = 0.3° com- 
pared to data from Ref. 0. The different curves are defined as in Fig. f|. Since 
it is impossible to use El at 1.4 GeV, and neither El nor E2 at 1.9 GeV, the cal- 
culations with constant amplitudes (dot-dashed line) are given at these energies. 
The calculations are reduced by the scale factors given in Sec. |3~6| . 
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Figure 6: Predictions of forward analyzing powers compared to the SPESIII 
measurements . Both the El (dashed line) and E2 (solid line) calculations are 
given. The experimental acceptance is incorporated in the theoretical results, 
(a) Average slope of vector analyzing power A! . The similarity in shape is more 
evident after shifting the El prediction by —2 rad -1 (dotted line). Data are 
obtained for a range of energies 1.116 < < 1.127 GeV, while calculations are 
done at 1.122 GeV. (b) Tensor analyzing power A yy at Td = 1.116 GeV. 
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